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We study the evolution of mixed scalar as well as spinor fields within the context of the classical field theory. 
The initial condition problem is solved and the fields distributions, exactly accounting for the initial conditions, 
are obtained for both scalar and spinor fields. In the system of two coupled fields we consider the special case of 
the initial conditions which are rapidly oscillating functions. It is demonstrated that the energy densities of the 
scalar fields and the intensities of the spinor fields coincide with the usual transition and survival probabilities 
of neutrino flavor oscillations in vacuum. 



1. INTRODUCTION 



The mixing between different flavor eigenstates of 
quarks as well as leptons was theoretically predicted 
more than fifty years ago (see Refs. [HQ). Since 
then a considerable progress has been achieved in both 
experimental and theoretical studying of elementary 
particles systems with mixing. One should mention 
the experimental discoveries of K and B mesons os- 
cillations (see, e.g., Refs. as well as the recent 
achievements in solar, atmospheric and reactor neu- 
trino experiments (see Refs. [!, H, @] ) • 

The theoretical description of neutrino flavor oscil- 
lations as the quantum mechanical evolution of two 
mixed flavor eigenstates was carried out in Ref. Q. 
More detailed quantum mechanical treatment of neu- 
trino flavor oscillations was given in Ref. @ where the 
wave packets approach was proposed. The quantum 
field theory was applied to the problem of neutrino 
flavor oscillations in Ref. in which the process of 
oscillations was explained as the propagation of neu- 
trino mass eigenstates described by internal lines of 
a Feynman diagram. The nonperturbative quantum 
field theory effects in neutrino flavor oscillations were 
discussed in Ref. [ll| where the authors considered 
the Fock spaces flavor and mass eigenstates and re- 
produced the Pontecorvo formula for the transition 
probability as well as obtained the corrections to this 
expression. We analyzed neutrino spin and flavor os- 
cillations on the basis of the classical field theory in 
Refs. [TH, [TH, It was shown that neutrino fla- 

vor oscillations in vacuum as well as in external fields 
can be described within the framework of the classi- 
cal field theory. Finally it is worth noticing that many 
other important references on the considered issue are 
presented in the review [l5[ devoted to the theory of 
neutrino flavor oscillations. 
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In this paper we continue to study neutrino flavor 
oscillations in vacuum within the context of the clas- 
sical field theory. Classical particles are represented 
with help of the first quantized fields in our approach. 
This terminology is borrowed from Ref. [l6j]. In Seed] 
we discuss the evolution of N coupled scalar fields. We 
solve the initial condition problem for this system and 
construct the fields distributions. The energy density 
of each field is calculated. It is shown that the ob- 
tained expressions are analogous to the formulae used 
in describing neutrino flavor oscillations in vacuum. 
In Sec. [3] the similar problem is solved for N cou- 
pled spinor fields. We demonstrate that the intensi- 
ties of the fermions are the same as the transition and 
survival probabilities of neutrino flavor oscillations in 
vacuum. We discuss our results in Sec. HJ 



2. EVOLUTION OF SCALAR PARTICLES 

Let us study the evolution of N coupled scalar fields. 
The Lagrangian for this system has the form 



A'* 



A'* 



(i) 



k=l 



i,k=l 
i^k 



where gik = g* ki are the coupling constants responsible 
for vacuum mixing, ip = (tpi, . . . , (Pn), and 



mt\vk\ 2 , 



(2) 



is the free field Lagrangian for the field ipk, rrifc is the 
mass corresponding to the field tpk- In Eq. |T]) we 
present the case of the complex fields (fk which can 
correspond to electrically charged particles. 

Following the results of our previous work [l3j we 
set the initial conditions problem to describe the evo- 
lution of our system {T]) and ([2]) . Supposing that the 
initial fields distributions 
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<Pi(r,0) = fi(r), <pi(T,0) = gi(r), 



(3) 
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are the given functions, one should find the fields 
distributions (p k (r, t) for any subsequent moments of 
time. It should be noted that it is necessary to set two 
initial conditions since the corresponding wave equa- 
tion - Klein-Gordon equation - is the second order 
differential equation. 

In order to solve the initial conditions problem we 
introduce the new set of the fields u k (r, t), which di- 
agonalize the Lagrangian (P) . They are usually called 
the mass eigenstates. The fields u k are related to the 
fields ifk by the matrix transformation, 



JV 



<Pi(r, t) = M ia u a (r, t). 



(4) 



a=l 



Note that the masses tr^ of the fields ipi are related 
to the masses m a of the fields u a by the following 
formula: 



TV 



(5) 



1 and gl). 

1311 we obtain the field 



which results from Eqs. ([T]), 
Using our previous work 
distribution of ipj (r, t) which is consistent with the 
initial conditions ([3]), 



N 



^(r,t) = £ M ja M~ l 

ia—1 

x J d 3 r'[j)„(r-r',f)/,(r') 
+ D a (r-r',t) 9i (r')], 



(6) 



where 



D a (r,t) 



_ f d 3 p ipr sin£ a t 



is the Pauli- Jordan function and £ a = y/p 2 + m 2 . 

It was demonstrated in Ref. jbt] that it was 
more convenient to use the momentum representation 
rather than the coordinate one. Therefore we rewrite 
Eq. ([6]) in the form 



N 



ia—1 



d 3 P 



(2tt) 3 
A»(p,*)$<(p)] 



J» T [Da(p,t)fi(p) 



(7) 



with help of the Fourier transforms of the Pauli- 
Jordan function, 

D a (p,t) = cos£ a t, D k (p,t) = SiU (8) 

Oct 



and the initial conditions, 

/i(p) = J d're-^Mr), 
ft (p) = / dhe-^giir). 



(9) 



It is interesting to demonstrate that the total energy 
is conserved in our system. Indeed, the total energy 
is given by the formula 



E(i) = J d 3 rT t o ° tal (r,i), 



where T®® ted is the time component of the energy- 
momentum tensor. Using the general expression for 



1 total' 



rpflV 

1 total 



TV 



dC 



d^ k + h.c. -Lif v , 



where the Lagrangian C is presented in Eq. ([T]), we 
receive the formula for E(t) in the form 

E(t)= J d 3 r(f](|^(r,t)| 2 + |V^.(r,t)| 2 (10) 
J ^ fc=i 



N 



i,k=l 



The initial value of the total energy can be obtained 
with help of Eqs. ©, © and dJUJ), 



E(0) = 



d 3 p 

J2^f 



N 



X>*(p)I 2 +p 2 i/*(p)I s 



fc=i 



N . 

m 2 fe |/ fc ( P )| 2 )+ ftfc/i (P)A(P) [• (11) 

i,fc=l ' 



i^k 



Using the explicit form for the field distribution 
(fj(r,t) [Eq. (J7J)] and the Fourier transforms of 
the Pauli-Jordan function [Eq. (jSJ)] one can rewrite 
Eq. (fTO)) in the following way: 



AT . ,3 

e(*)= x; ^M-^y — p 



ika— 1 

X 



(2tt) 3 

[^(p)5*(p) + ^/i t (P)A(p) 



(12) 



In deriving Eq. (|12[) we use the unitarity of the matrix 
Mj a , i.e. (Af _1 ) a i = M* a , and the property of the 
masses eigenstates u a , 



N 



m 2 a 6ab = J2 m l M kaM k b + 9*kM* a M kb . 



k=l 



i,k=l 
i^k 
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With help of the identity 



N 



m; 



fife ■ 



which generalizes Eq. ([5]) we can transform Eq. 
to the form which coincides with Eq. (fTTj) . i.e. we 
demonstrate that E(i) = E(0). This explicit calcu- 
lation shows that the total energy in the system of 
N coupled scalar fields is conserved independently of 
initial conditions. 

It is very difficult to analyze Eq. (O for arbi- 
trary functions fi(p) and 5i(p). Thus we pick out 
the special form of the initial conditions (see also 
Refs. [13, El), 



It is also necessary to fix the "amplitudes" A4. They 
can be chosen in following form, A\ — and A 2 = 
1 . This choice of the " amplitudes" has rather simple 
physical meaning: the first field is absent initially and 
the second one has the plane wave field distribution. 

Using Eqs. Q, (JUJ) and (TT7D we g et the neld dis- 
tributions of ipi t 2\ 



<pi(r,t) 



V2M) 



cos sin 9e u 



x (cos[Si(u))t\ — cos{£ 2 {uj)t]), 



(18) 



x (sin 2 9 cos[£t(w)t] + cos 2 9 cos[£ 2 (oj)t\), 



Mr) = 



9i(r) = 0, 



(13) 



w here Aj i s the "amplitude" of the function /j, = 
\J u)" 2 + m 2 and u> are the initial energy and momen- 
tum of the field (pi. Note that is the normal- 
ization factor; its value will be clarified below [see 
Eq. (p~6|) ] . Using Eq. ([9]) we obtain the Fourier trans- 
forms of the initial conditions, 



where E a ((jS) = \Juj 2 + m\. To analyze the obtained 
expressions we consider the limiting case of high initial 
frequency, oj ^> mi,2, which corresponds to ultrarela- 
tivistic particles. 

With help of Eqs. ([T5]) and (JTHJ) as well as in the high 
frequency approximation we obtain the expressions for 
the energy densities of the fields <p\ 2 in the following 
form: 



f i (p) = ^L(2*) 3 5 !i (u,-p), 9i (p) = 0. (14) 



The physical quantity measured in an experiment 
is not a field distribution. It was demonstrated above 
that the total energy of the system is conserved. 
Therefore we assume that one can detect the en- 
ergy density of a single scalar field. Constructing the 
energy-momentum tensor for the single field ipi and 
taking its time component we obtain for the energy 
density, 



«,(r,i) =T°°[ ¥ > i (r,t)] 

=|<^| 2 + \v<pi 



m 2 |^| 2 - 



(15) 



Note that here we assume that there would be no mix- 
ing between different fields ipi. It can be verified by 
means of direct calculations that the following expres- 
sion: 



J5i(r,0) = \A t \ 2 ^ 



(16) 



results from Eqs. (|13p and (fTS"]) . Therefore the nor- 
malization factor in Eq. (TIB"]) was chosen to get the 
correct value of the initial energy density. 

Without losing generality we can discuss now the 
evolution of two scalar fields. In this simple case the 
matrix Mj a can be parameterized with help of only 



one mixing angle 9, 



Ma 



( cos 9 
{ sin 9 



— sm9 
cos 9 



(17) 



=w <^ sin 2 (26>) sin 2 [A(w)t] + O 



> (19) 



fait) =u)\l- sin 2 (26») sin 2 [A(w)£] + O 



where 



Am 2 ^ / m 2 

-n— + °\-§- 



Here we use the common notation Am 2 = m\ — m\. 

Now we can introduce the probability which corre- 
sponds to the transitions from the second eigenstate 
tp2 to the eigenstate ifk with k — 1,2: P^ft CO = 
9)k(t)/u>. Using Eqs. ([IT))) we obtain for the probabil- 
ities, 



^(scalar) 



(t) = sin 2 (26 1 ) sin 2 



Am 

4w 



(20) 



^(scalar) (t) =J _ ^ ^ A (21) 



The function P 2 ( !^ lar) (i) in Eq. (2DJ) is usually 
called the transition probability whereas the function 
P 2 (^ lar ) m Eq pjj) - the survival probability. It 
should be noted that Eqs. (|20p and (j2Tj) reproduce the 
common quantum mechanical expressions for transi- 
tion and survival probabilities of neutrino flavor oscil- 
lations in vacuum. 
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3. EVOLUTION OF SPINOR PARTICLES 



where 



Now let us discuss the evolution of N coupled spinor 
fields. The Lagrangian for this system has the form 



JV 



TV 



= E £oK) - E 9tk 



ViVk-, 



k=\ 



i,k=l 
i^k 



where v = (i>\,..., vn) and 

^o(ffc) = vu^dn - m k )v k , 

is the free field Lagrangian for v k and rrifc are the 
masses of v k . 

Analogously to Sec. [2] (see also Refs. [H, we 
set the initial condition problem for the considered 
system, 



i/i(r,0)=&(r), 



(22) 



where £i(r) are the given functions. It should be men- 
tioned that we have to set only one initial condition 
in Eq. (f2"2"|) since the Dirac equation is the first order 
differential equation. Now it is necessary to find the 
fields distributions for ^(r, i) at t > which would be 
consistent with Eq. (l22l) . 

With help of Refs. |l3l.[l4j we obtain for Vj(r, t) the 
following expression: 

N 



v j {r,t) = Y,M ja M-? 
=i 

J dVS a (r'-r,t)(-i 7 )&(r'), (23) 



ia—l 
X 



where 

S a (r,t) = [Yfd M +m a )D a (r,t), x" = (t,r), 

is the Pauli- Jordan function for a spinor field and m a 
are the masses of mass eigenstates which are related 
to rrii by the formula, 



N 



m, 



a=l 



It should be noted that Eq. |23J is in agreement with 
the initial conditions in Eq. (|22)). 

As it was demonstrated in Sec. [5] we have to rewrite 
Eq. (|23[) in the momentum representation in order to 
avoid computational difficulties. Finally we get 

jV 



i (r,t) = Y,M ja M- i l 



d 3 p 

(2tt)3 



e^S a (-p,t)(-h°Mp), (24) 



S a (—P, t) = COsSat 

.s'm£ a t 

i— p — (ap + (3m a ) 



(25) 



and 



are the Fourier transforms of the Pauli-Jordan func- 
tion and the initial conditions. Here we adopt the 
common notation for the Dirac matrices a = 7 ° 7 and 
/3 = 7° 

It is interesting to verify that in the classical field 
theory approach the total intensity of all the fields v k 
is conserved. We assume that 



JV 



(())=£ /d 3 r|6(r) 



i=l 
JV 

E 



d 3 P 



I — 1 v y 



I6:(P)| 2 , 



(26) 



has the given value. Now we are interested in the 
following quantity: 



N 

J(t)=E _/d 3 rh(r,f)p 



With help of Eq. (|24p this expression can be rewritten 
in the form 

TV 

3(t)= ]T [M^M-^M^M-,} 

jikab—1 

. , d3 q d 3 re i( q -P)r 



(2tt) 3 (2tt) 3 

xCl(p)[5 o (-p,f)(-i 7 )] t 

x ,S 6 (-q,t)(-i 7 %(q). (27) 

Again using the fact that the matrix Mj a is the uni- 
tary one, we obtain for Eq. (f2"7| , 

N 
ika—1 



(2tt) 3 ?l 



x I ^e!(p)[^(-P,0(-i7°)] t 
xS a (-p,t)(-i 7 )&(p). (28) 



With help of Eq. (125[) one can show that 

[^(- P ,i)(-i 7 )]t5 a (-p,t)(-i 7 °) = l. 
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Therefore, taking into account Eq. (|26|) . we represent 
Eq. (J2HD in the form 



N 



d 3 P 

(2tt) 3 



|C4p)| 2 =3(0), 



(29) 



which shows that the the total intensity is conserved. 

In the following we consider the evolution of only 
two coupled spinor fields to make the results more il- 
lustrative. Thus we assume that the initial conditions 
are: £i(r) = and £2(1") = e la ""£o. The discussion of 
such a choice of initial conditions is given in Sec. El 
The mixing matrix Mj a can be also represented in the 
same form as in Eq. (JT7J) . Using general Eq. (|24p one 
obtains for the fields distributions of v\ 2 ; 



Eqs. ([31]) we express the probabilities in the following 
way: 



^(spinor) ( ^ =sm 2 (2 )sin 2 



Aw 



^(spinor) 



(t) =1- sin 2 (26») sin 2 



Am 2 

4uj 



(32) 
(33) 



One can observe that these expressions for transi- 
tion (|3"2"]1 and survival (|33[) probabilities coincide with 
the analogous Eqs. (20) and (21) for the coupled scalar 
fields. Note that Eqs. ([3"2"| and ([3"3")) are consistent with 
Eq. (|2"5)) . The obtained probabilities formulas are the 
same as in the quantum mechanical treatment of neu- 
trino flavor oscillations in vacuum. 



vi(r,t) =e la,r sin cos 0{ cos[£i(u>)t] — cos[£2(w)i] 
. sin[f 



(aw + f}m±) 
(au + f3m 2 ) !>£o, 



.sin [£2 (<*>)*] 

1 



u 2 (r,t) =e l " r \ sin 2 0cos[£i(w)t] + cos 2 cos[£ 2 (u)t) 



1 sm 2 — i , / J (au + /3mi] 



i cos 6 



2 sm[£ 2 (u)t] 



£ 2 (w) 



(aw + Mzj ko. (30) 



The most interesting case is the high frequency ap- 
proximation of the initial conditions given in Eq. (|2"2"J) , 
i.e. uj mi^- Therefore, Eqs. ([3TJ]) can be rewritten 
in the form 

v x ( r , t) = - e [ujr sin 26> sin[A(w)t] 

x {sin[<r(a;)i] + i(an) cos[cr(a;)rj]} £0 
+ 0(m a /u), 

v 2 (r,t) =e iwr {sin[A(w)*] cos 20 - i(an) cos[A(w)t]} 
x {sin[a-(ti>)t] + i(an) cos[cr(w)t]} ^ 
+ O(m„/o;), (31) 

where 



a(w) = + 



We showed above that the total intensity of all the 
fields z/fc, which is proportional to a squared fields 
distributions, is conserved in the considered system. 
Therefore we suppose that the measurable quantity 
of a classical spinor particle is the intensity, Ik(t) = 
(r, t)\ 2 . Analogously to Sec. Owe introduce the 
probabilities for transitions from the second eigenstate 
to the eigenstate k = 1, 2 as P, 2!^^°^ (t) = Ik(t)- Using 



4. CONCLUSION 

In conclusion we note that the evolution of coupled 
scalar as well as spinor fields has been described in 
frames of classical field theory. First we have con- 
sidered the system of N arbitrary coupled complex 
scalar fields and solved the initial conditions problem 
for this system. The general fields distributions con- 
sistent with the initial conditions have been found. 
Note that these expressions exactly took into ac- 
count the Lorentz invariance and were valid in (3+1)- 
dimensional space-time. The energy conservation law 
has been analyzed for the considered system. Then we 
have discussed the simple case of two coupled fields 
and have constructed the energy densities of each 
scalar field. It has been demonstrated that in ultrarel- 
ativistic limit these expressions were analogous to the 
transition and the survival probabilities of neutrino 
flavor oscillations in vacuum. The calculations in the 
presentpaper have refined the results of our previous 
work [131 ] where analogous problem was analyzed in 
(l+l)-dimensional space-time. 

The case of N coupled classical spinor fields has 
been also considered in the present work. We have 
formulated the initial condition problem and con- 
structed the general expressions for the fields distri- 
butions which were valid in (3+l)-dimensional space- 
time and were Lorentz invariant. The total intensity 
of all coupled spinor fields has been computed. It has 
been demonstrated that this quantity was conserved 
independently of the initial conditions. Then we have 
discussed the case of rapidly oscillating initial condi- 
tions which corresponded to ultrarelativistic particles. 
In this case as well as for only two coupled fermions we 
have obtained the intensities of the fields in question. 
The expressions for the intensities were also revealed 
to coincide with transition and survival probabilities 
of neutrino flavor oscillations in vacuum. It should be 
mentioned that transition and survival probabilities 
for coupled classical spinor fields turned out to be the 
same as the analogous formulae for scalar fields. 
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